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ABSTRACT 

This  study  considers  the  problem  of  designing  a  flexible  multi-line  facility  that  manufact  ures 
a  variety  of  products  in  large  volumes  with  stable  demand  rates.  These  products  have 
similar  processing  requirements  in  that  they  visit  the  various  manufacturing  stages  in  the 
same  sequence.  Each  stage  on  a  given  line  performs  a  set  of  predetermined  tasks,  and  it 
comprises  multiple  identical  CNC  machines.  While  these  tasks  require  similar  processing 
capabilit  ies,  the  actual  tasks  done  and  their  processing  times  are  product-specific.  This 
paper  considers  the  special  case  in  which  there  is  only  one  stage  of  manufacture  -  this  case 
arises  as  an  important  subproblem  while  solving  the  more  general  multi-stage  problem. 

We  consider  two  versions  of  the  problem,  one  in  which  the  demand  of  any  product  can 
be  split  across  multiple  lines,  and  the  other,  restricted  version  in  which  each  product  is 
assigned  to  exactly  one  line.  For  the  demand  splitting  case,  we  discuss  the  properties  of 
the  optimal  solution,  and  present  a  polynomial  time  algorithm  for  solving  it.  The  latter 
problem  is  shown  to  be  NP-hard  in  the  strong  sense.  We  discuss  several  relaxations  of  this 
problem,  and  develop  lower  bounds  based  on  these  relaxations.  We  also  construct  an  implicit 
enumeration  based  solution  algorithm  for  this  problem. 


1      Introduction 

This  study  considers  the  flexible  multi-line  design  (FMD)  problem  in  a  serial,  multi-stage 
manufacturing  facility.  This  facility  processes  a  variety  of  products  in  large  volumes  with 
stable  demand  rates.  These  products  have  similar  processing  requirements  in  that  they  visit 
the  various  manufacturing  stages  in  the  same  sequence.  Each  stage  on  a  given  line  performs 
a  set  of  predetermined  tasks,  and  it  comprises  multiple  identical  CNC  machines.  While  these 
tasks  require  similar  processing  capabilities,  the  actual  tasks  done  and  their  processing  times 
arc  product-specific.  The  flexible  CNC  machines  can  switch  from  one  product  to  another 
with  negligible  changeover  time.  Products  are  moved  from  one  machine  to  another  in  small 
cassettes;  the  size  of  these  transfer  batches  is  quite  small. 

Each  line  is  paced.  The  cycle  time  for  a  given  line  is  determined  by  the  products,  in  particular 
the  bottleneck,  hereafter  the  pivot  product,  assigned  to  it.  While  there  may  be  economic 
incentives  in  having  multiple  lines  in  order  to  reduce  the  idle  times  at  individual  stages,  the 
benefits  of  doing  so  need  to  be  traded  off  against  the  fixed  cost  of  providing  the  line.  The 
decision  to  be  made  at  the  design  stage  is  to  determine  the  number  of  lines  to  be  formed  and 
to  find  the  number  of  machines  to  be  assigned  to  each  stage  on  each  line  such  that  the  total 
investment  in  lines  and  workcenters  required  to  meet  the  given  demand  rates  is  minimized. 

This  problem  is  motivated  by  the  manufacturing  facility  of  one  of  the  major  auto  companies 
that  produces  fuel-supply  systems.  This  problem  also  arises  in  printed  circuit  board  manu- 
facture (Farber  et  al.  1988).  The  FMD  problem  generalizes  the  mixed-model  line  balancing 
problem  considered  in  Miltenberg  (1989),  Inman  and  Bulfin  (1991),  Kubiak  and  Sethi  (1991). 
Miltenberg  and  Sinnamon  (1989),  and  Okamura  and  Yamashita  (1979)  to  multiple  lines  and 
parallel  workcenters  available  at  a  given  stage.  In  this  paper,  we  consider  the  special  case 
in  which  there  is  only  one  stage.  This  special  case  arises  as  an  important  subproblem  while 
solving  the  general  multi-stage  problem. 

I  he  paper  is  organized  as  follows.  We  formulate  the  FMD  problem  in  §2.  We  consider  two 
versions  of  the  problem,  one  in  which  the  demand  of  any  product  can  be  split  across  multiple 
lines,  and  the  other,  restricted  version  in  which  each  product  is  assigned  to  exactly  one  line. 
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We  show  that  the  latter  problem  is  NP-hard  in  the  strong  sense.  §3  deals  with  the  demand 
splitting  case.  We  discuss  the  properties  of  the  optimal  solution,  and  present  a  polynomial 
time  algorithm  for  solving  it.  In  §4,  we  discuss  the  restricted  problem,  and  present  several 
lower  bounds  for  it  in  §5.  An  implicit  enumeration  based  algorithm  for  solving  this  problem 
is  given  in  §6.  We  conclude  in  §7  with  a  summary  of  the  main  results  of  this  paper. 

2      Description  of  the  Problem 

First  we  give  the  notation  used  in  the  paper. 

Af  :  the  set  of  products,  and  \Af\  =  N 

F\  :  fixed  cost  per  line 

F2  :  fixed  cost  per  machine 

Pj  :  processing  time  of  product  j,    j  £  J\f 

dj  :  daily  demand  of  product  j,    j  €  Af 

A  :  available  time  per  day  on  any  machine 

B  :  a  large  number 

n3  :  number  of  machines  required  on  a  line  with  pivot  j 

Xji  =  the  fraction  of  product  i  assigned  to  the  line  with  pivot  j 

1,    if  a  line  is  opened  with  pivot  j 

Vj    =     \ 

0,    otherwise. 

2.1.     Problem  Formulation 

Throughout  this  paper  we  assume  that  the  products  are  numbered  such  that  if  i  <  j,  i,  j  £ 
TV,  then  pi  >  pj.  Also,  let  Ji  =  {j\p3  >  pt}.  The  product  assignment  problem  can  now  be 
stated  as: 


PI 


P 
Minimize      Zx  =  Yl^Vi  +  F^nj)  t1) 


subject  to 

£*0-  =  l,     i€tf  (2) 

j€.7, 

Pife^iJ  <**i,     i€^  (3) 

n3<Byr     je.V  (4) 

xo€{0.1},     i.jGA'  (5) 

I) j  e  {0,  l};n;  >  0,    integer,     ;  G  A'  (6) 

Equation  (2)  insures  that  a  product  is  assigned  to  exactly  one  line.  Constraint  (3)  specifies 
that  the  cycle  time  for  all  products  assigned  to  a  line  with  pivot  j  is  pj.  and  the  number 
of  required  machines  is  computed  accordingly.  Constraint  (4)  insures  that  the  fixed  cost  of 
opening  a  line  is  accounted  for.  Constraints  (5)  and  (6)  specify  the  nature  of  the  variables. 
The  demand  splitting  version  of  the  product  assignment  problem  is 

P2 

P 

Minimize      Z2  =  ^2(Fiyj  +  ^2">)  (") 

;=i 

subject  to 

(7) -(4),  and 
0<xti<l,      ij€Af  (8) 

Theorem  1.    Pi  is  NP-hard  in  the  strong  sense. 

PROOF:  PI  is  clearly  in  NP.  We  show  that  the  3-PARTITION  problem  which  is  known 
to  the  NP-hard  in  the  strong  sense  (Garey  and  Johnson  1979)  is  reducible  to  Pi.  Consider 
an  arbitrary  instance  of  3-PARTITION  given  by  a  set  Q  of  3<?  elements  of  size  st  for  each 
i  6  calQ.  and  a  positive  integer  B  such  that  i)  Bj\  <  s,  <  B/2,  Vz  £  Q.  and  ii)  2I.ee5'  = 
qB.  The  recognition  problem  is  stated  as:  Is  there  a  partition  of  Q  into  q  disjoint  subsets 
1  Qi,  Qi Qm)  such  that  E.60;  -s>  =  B^  for  *  ^  J  ^  <?? 


The  equivalent  instance  of  PI  is  constructed  as  follows: 

\Af\     =  Aq 

Pj    =  h  3  =  1,2,..., 9 

=  1,  j  =  q  +  1,^  +  2,..., iq 

dj     =  qB,  j  =  1,2,  ...,q 

=  Sj-q,    j  =  q  +  1,9  +  2,...,  4q 

A     =  qB  +  B 

F1    =  0 

Given  an  instance  of  3-PARTITION,  this  instance  of  Pi  can  be  constructed  in  polynomial 
time.  We  show  that  Zl  <  F2q(q  +  l)/2,  if  and  only  if  3-PARTITION  has  a  solution.  For 
ease  of  exposition,  denote  the  set  of  products  j  =  1, 2, . . . ,  q  by  Af\,  and  let  V2  =  N  —  N\- 
First  suppose  that  3-PARTITION  has  a  solution.  The  required  solution  for  PI  is  constructed 
by  forming  q  lines  such  that  each  product  j  £  A/*i  is  a  pivot.  Note  that  the  number  of  machines 
required  individually  for  any  product  assigned  to  a  line  with  pivot  j  is 

jqB 


N     = 
and  the  resulting  idle  capacity  is 


(qB  +  B) 


=  J  (9) 


j{{qB  +  B)  -  qB}/j  =  B  units.  (10) 

The  total  number  of  machines  required  under  this  assignment  for  products  in  M\  is  Yll=i  I  — 
q(q  +  l)/2.  Clearly,  if  3-PARTITION  has  a  solution,  M2  can  be  partitioned  into  q  disjoint 
subsets  7^21,  A/22,  •••  ,A/2g  such  that  YL]£V2i  dj  —  B,  I  =  1,2, . . .  ,q.  Each  of  the  q  lines  is 
assigned  one  of  these  subsets  to  give  the  desired  solution. 

Next,  we  show  that  if  z  <  F2q{q  +  l)/2  under  any  given  assignment  a  (say),  then  3- 
PARTITION  has  a  solution.  First,  note  that  because  p,  >  pj,  i  gMiJ  £  A/2,  there  is  at  least 
one  line  in  a  that  has  a  pivot  which  belongs  to  M\.  Suppose  that  this  solution  has  L  +  T  lines, 
of  which  the  pivots  for  the  first  L  lines  are  ji,J2,  •  •  •  ,Jl,  such  that  ji  £  V\,l  —  1,2, . . . ,  L. 
Clearly,  L  <  q.  Let  Mi  be  the  number  of  machines  in  line  /,    /  =  1,  2, . . . ,  L,  and  let 

Ci     =     {j\j  £  Af\,j  is  assigned  to  /} 


7/     =     |£/|- 

First  consider  the  assignment  of  products  in  Pi  to  these  L  lines. 

Lemma  1.    £/={/}.    I  =  1,2 L,  and  hence,  L  =  q. 

PROOF:  Let  il5 i2, . . . ,  <->,  (=  /)  be  the  products  assigned  to  any  given  line  /  under  a.  Consider 
an  alternative  assignment  a'  in  which  products  *i , «2>  •  -  * » *-vi— l  are  pivots  for  their  individual 
lines.  From  (1),  the  total  number  of  machines  required  under  this  configuration  for  these  7/ 
products  is 


<     (7/ -!)(/-!)  +  / 

Note  that  (3)  is  satisfied  as  an  equality  only  if  either  C\  =  {/},  or  Ci  =  {I  —  1.  /}.  Now 

IqB-n 


(ID 


Mi    = 


qB  +  B 

,        hi 
hi  - 


>  h/i  —  (7/  —  1 ),    because  /  <  q  +  1 

>  /+(7/-l)(/-l) 


Hence, 


L  L 

£  A/,  >  £  M;  =  1  +  2  +  . . .  +  q  =  q{q  +  l)/2 


Given  that  Zi  <  F2<?(r/  +  l)/2,  it  follows  that 


£Af,  =  «(g+l)/2 


and  £/  =  {/},  or  £/  =  {/  —  1 ,  /}.  V/.  Now  suppose  that  for  any  line  /,  Ci  =  {I  —  1.  /}.  The 
number  of  machines  required  for  this  line  is  (2/  —  1),  and  the  capacity  available  for  products 
in  Mi  is 

CAP(l)  =  [{21  -  \)(qB  +  B)  -  2lqB]/l  <  IB  units.  ( 12) 


From  (10)  and  (12),  the  total  capacity  available  on  all  L  lines  for  products  in  ^f2 

L 

J2  CAP{1)  <  2B  +  {L-  2)B  <  qB. 
1=1 

But  the  minimum  capacity  required  for  products  in  Af2  =  J2%q+i  dj-q  —  Xlf=i  si  —  qB.  Hence, 
for  the  only  feasible  configuration  under  a  is  £/  =  {/},  /  =  1, 2, . . . , L,  and  L  =  q. 
Now  consider  the  products  in  A/2.  Given  that  there  are  q  lines,  each  with  idle  capac- 
ity B,  a  feasible  solution  can  exist  only  if  A/2  can  be  partitioned  into  q  disjoint  subsets 
A/"2i,A/,22,  •  • .  1A/'2q  such  that  Zj€V2l  dj=B,  I  =  1,2,. . .  ,g.  This  implies  that  3-PARTITION 
has  a  solution.  □ 


3      Demand  Splitting 

We  first  address  the  problem  in  which  the  demand  of  a  product  can  split  across  multiple 
lines.  Consider  an  arbitrary  solution  a  with  L  lines.  Let  V  =  {j\j  E  Af,yj  =  1}  denote 
the  set  of  pivots,  and  C3,j  G  V  denote  the  set  of  products  assigned  to  pivot  j .  To  simplify 
exposition,  we  use  (/)  to  denote  the  pivot  for  line  /,  /  =  1, . . . ,  L,  and  [k],  k  E  V  to  denote 
the  line  for  which  k  is  the  pivot.  The  number  of  machines  required  on  line  /  is  given  by 

(diXn)p{i) 


Mi  = 


A 


and  Z2(a)  =  FXL  +  F2Eti  Mi. 


Proposition  1.    If  F\  >  0,  then  in  an  optimal  solution,  p^  >  p(/j  if  k  >  I,  k,l  G  V . 

PROOF:  From  the  product  numbering  scheme,  p^  >  pay  Suppose  that  a  is  an  optimal 
solution  to  P2  that  does  not  have  the  above  property.  Then  p^)  —  P(i)-  Construct  an 
alternative  solution  a'  in  which  line  [k]  is  merged  with  line  [/],  and  the  assignments  on  other 
lines  are  the  same  as  in  a.  Now 


Z2(a)  -  Z2(a')    =    2Fl  +  F2< 

-Fx  -  F2 
>    0 


P(l)  E«G£,  dtXU 


A 


+ 


P{k)  J2teck  dtxtk 
A 


P(i) 


[Y,tec,  dtXu  +  Yteck  dtxtk) 


Ix-caiisc  of  t he  known  inequalil y 

[a +  61  <  \a]  +  \b]  (13) 

and  the  fact  that  F\  >  0.  and  /?(*.)  =  p^/j.  Hence,  a  cannot  be  optimal,  and  the  proof  is 
complete.  □ 

It  is  easy  to  see  that 

Corollary  1.    //  F\  =  0,  then  it  is  optimal  to  have  a  single  line  with  product  1  as  the  pivot. 

Henceforth,  we  restrict  ourselves  to  the  cases  in  which  F\  >  0.  Throughout  this  paper,  we 
assume  that  the  lines  are  numbered  in  the  decreasing  order  of  their  cycle  times. 

Proposition  2.    (Sequential  Assignment  Property):  There  exists  an  optimal  solution  to 
Pi  in  which,  for  any  non-pivot  product  i,  xul  G  {0, 1}  for  u  €  P,  u  <  i.   Moreover,  if 
Xji  =  1.  then  Xjq  =  1  for  q  —  j  -f  1,  j  +  2, . . . ,  i '■  —  1,  i  as  well. 

PROOF:  Let  a  be  an  optimal  solution  to  P2  that  does  not  satisfy  the  sequential  assignment 
property.  Let  Z\  denote  the  optimal  solution  value,  and  V  denote  the  set  of  pivots  in  0\. 
Let  t  be  a  non-pivot  product  that  is  assigned  to  J,  L  >  J  >  1  lines  in  a .  Let  Vt  denote  the 
set  of  pivots  corresponding  to  these  lines.  Without  loss  of  generality,  we  assume  that  these 
are  the  first  J  pivots  in  "P,  i.e.,  xu  >  0.  for  /  =  ( 1 ),  (2), . . . ,  ( J).  Note  that 

Pt<P(i),    vievt  (14) 

First,  we  establish  the  following  result. 

Lemma  2.     There  exists  an  alternative  optimal  solution  a'  with  the  property  that  eitht  r  i) 
x{j)t  =  Yiq=i  x(q)t'  or  H)  t  *5  a  pivot  for  line  J,  where  the  "prime"  denotes  the  assignment 
variables  in  a' . 

PROOF:  The  proof  is  by  induction.  We  show  that  if  the  result  is  true  for  the  first  k  lines. 
then  it  is  also  true  for  the  first  k  +  1  lines;  first,  we  show  that  the  result  holds  for  k  =  2  lines. 


Let  D\{D2)  denote  the  total  number  of  units  allocated  to  line  1  (2).  Construct  a'  from  a  as 
follows:  put  X/jw  =  0,  and  x',2u  =  £(i)t  +  X(2)t-  This  modification  shifts  6t  =  dtx^t  units  of 
demand  from  line  1  to  line  2.  Replace  these  units  on  line  1  by  products  currently  assigned 
to  line  2  considered  in  the  increasing  order  of  their  index.  At  the  end  of  this  step, 

PiD1 


Mi(<t')  =  M1(<t)  = 


A 


From  the  construction  of  cr',  it  is  clear  that  if  d(2)X(2)(2)  >  $t,  then  (2)  continues  to  remain  a 
pivot  for  line  2,  and 

M2(a')  =  M2(a)  = 


P(2)A 


A 

otherwise,  the  new  pivot  for  line  2  is  j  such  that  p3  <  p(2),  and 

-pjD2 


M2(a')  = 


A 


<  M2{a). 


In  either  case,  Z2(a')  <  Z2(a)  because  the  number  of  lines  is  the  same  in  both  a  and  a'.  If 
j  <  t,  then  X/2w  =  X(i)t  +  x^)f  Otherwise,  j  =  t  is  the  pivot  for  line  2.  Note  that  j  >  t  is 
not  possible  because  this  implies  that  x<2)t  =  0  which  contradicts  our  assumption. 

Suppose  that  the  result  holds  for  k  =  5,  i.e.,  we  have  an  alternative  optimal  solution  a'  such 
that  either 

Case  1:  x'(s)t  =  £,aq=lx(q)t,  or 
Case  2:  t  is  a  pivot  for  line  s. 

For  Case  1,  we  next  consider  lines  s  and  5  +  1;  the  above  arguments  can  be  used  to  show  that 
the  result  holds  for  k  =  s  -f  1  also.  For  Case  2,  consider  lines  5  and  s  +  2.  The  arguments 
used  for  k  =  2  can  be  used  to  show  that  either 

Case  3:  x[s+2)t  =  £J=1  x(q)t  +  x(,+2)t,  or 
Case  4:  £  is  a  pivot  for  line  s  +  2. 

However,  Case  4  implies  that  /  is  a  pivot  for  lines  5  +  1  and  s  +  2.  From  proposition  3, 
it  follows  that  cr  is  not  optimal  for  P2  which  contradicts  our  assumption.  Hence,  if  Case 
2  holds  for  k  =  5,  then  Case  3  must  be  true  for  k  =  s  +  1.  Renumbering  the  lines  in  the 
decreasing  order  of  their  cycle  times  insures  that  t  is  the  pivot  for  line  5  +  1.  □. 
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Applying  Lemma  2  to  each  non  pivol  product  thai  is  assigned  to  multiple  lines  in  g  eventually 
yields  an  alternative  optimal  solution  n'  thai  satisfies  the  first  part  ol  proposition.  If  n'  does 

not  satisfy  the  second  part  of  the  proposition,  then  there  has  to  be  at  least  on<-  producl 
t,j  <  t  <  i  —  1  that  is  produced  on  line  [k],  k  ^  j.  Clearly,  k  <  t  <  i..  Let  S  =  <l,.r.: 
be  the  demand  of  t  allocated  to  line  [k\.  If  k  <  j,  then  construct  a"  from  <i'  by  shifting  i 
units  of  t  to  line  [j],  and  replace  these  units  on  line  [k]  by  products  currently  assigned  to 
line  [j]  considered  in  the  increasing  order  of  their  index  starting  with  j .  If  S  <  djXjj,  then 
j  continues  to  remain  a  pivot,  otherwise  its  entire  demand  is  absorbed  by  line  [k]  and  is 
replace  by  a  pivot  q\q  >  j.  with  p?  <  pj.  In  either  case,  Z2(cr")  <  Z2(t/),  and  therefore,  a" 
is  optimal. 

[f  fc  >  j,  then  construct  cr"  by  shifting  6  =  Ylq=t+\  CK  units  of  demand  corresponding  to 
products  t  +  1  through  i  from  line  [j]  to  line  [&],  and  replace  these  units  on  line  [j]  with 
products  currently  assigned  to  line  [k]  considered  in  the  increasing  order  of  their  index 
starting  with  k.  As  before,  it  follows  that  Z2(<t")  <  Z2(cr'),  and  therefore,  a"  is  optimal. 
Repeating  these  steps  whenever  required  yields  the  solution  a"  that  is  optimal  to  P2  and 
that  satisfies  the  condition  stated  in  proposition  3.  □ 

An  immediate  consequence  of  the  sequential  assignment  property  is  that  if  j  is  a  pivot  in  an 
optimal  solution,  then  it  is  assigned  to  at  most  two  lines,  namely  [j '  —  1]  and  [j].  Now  consider 
line  1.  Clearly  product  1  must  be  a  pivot  for  this  line.  Consider  an  arbitrary  solution  in 
which  .rlu  =  1  for  u  =  1 i,.  The  number  of  machines  required  at  line  1  is 


Mu  = 


A 


The  capacity  remaining,  hereafter  the  remnant,  at  line  1  after  this  assignment  is 


r^,  = 


-5X 


Pi  u=l 

Clearly,  it  is  optimal  to  use  this  capacity  for  producing  product  i  ■+■  1,  such  that  .rll  +  1  >  0. 
In  the  general  case,  the  remnant  available  at  line  /  will  be  used  for  producing  (/  +  1). 
Consequently,  the  number  of  machines  required  by  a  line  with  pivot  j.  and  rJU  =  1  for 
u  —  j '  +  1 k  is 

Pj  (eIu.,  d*  -  m)-u-\) 

Mjk 


and 


'•iit  = 


AMj, 
Pj 


£  du  ~  r(j)-u-i     ■ 


<u—j 


Note  that  rjk  <  A/pj,  hence,  it  is  strictly  less  than  one  machine's  capacity  on  line  (j).  Also, 
because  p3  >  P(j)+i,  it  is  less  than  one  machine's  capacity  on  subsequent  lines  as  well.  The 
marginal  increase  in  the  number  of  machines  required  on  line  (j)  for  producing  the  k  +  \th 
product  is 


mjik+i  =  Mj,k+i  -  Mjk  = 


Pj  (4+i  -  rjk) 
A 


(15) 


P2  can  be  posed  as  the  shortest  path  problem  on  graph  G  =  (V,«4)  shown  in  Figure  1.  In  G, 
node  Vij  represents  an  assignment  in  which  product  j  is  produced  on  line  with  pivot  i.  Note 
that  node  Vij  is  feasible  only  if  pj  <  pt;  hence,  the  upper  triangular  nature  of  this  graph.  We 
append  a  dummy  node  T  at  the  end  to  represent  an  artificial  product.  The  optimal  solution 
to  P2  corresponds  to  the  shortest  path  in  G'  from  Vn  to  T.  Let  A"?  denote  the  arc  leading 
from  Vij  to  Vuv.  A  can  be  partitioned  into  subsets  7i,  B  and  T  such  that  7i  U  B  U  T  =  A, 
and  Tif)  B  =  Hd  T  =  BV\T  =  0.  7i  is  the  set  of  horizontal  arcs  (h-arcs)  of  the  form  A\-+  , 
while  B  is  the  set  of  backward  arcs  (b-arcs)  of  the  form  A??  where  u  <  i.  T  comprises  the 
forward  arcs  (f-arcs).  Without  any  loss  of  generality,  we  assume  that  all  arcs  leading  to  T 
are  f-arcs. 

INSERT  FIGURE  1  HERE 

From  proposition  3  it  follows  that  any  path  that  includes  a  b-arc  is  not  dominant.  Further- 
more, the  only  f-arcs  that  need  to  be  considered  are  of  the  form  A\j  ,J+  ,  i.e.,  those  arcs  that 
are  incident  on  a  pivot.  The  cost  of  such  an  arc  is 

Pj+i  (^i+i  —  rij) 


eg1'*1  =  Fl  +  F2 


(16) 


and  the  cost  of  the  h-arc  A)j       is 


.«d+i  _ 


*j 


=  F2 


Pi  {dj+i-  rtJ) 
A 


(17) 


Note  that  there  are  ;'  f-arcs,  of  the  form  AJ+1,i+1,  AJ+lj'+1, . . .  ,AJ+lj'+1  that  are  incident 
upon  Vj+ij+i.   From  16,  it  follows  that  the  costs  of  these  arcs  depend  upon  ry,  and  hence, 
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they  are  different.  However,  because  rtJ  <  A/pi  <  Pj+\,  for  i  =  1,. . .  ,j, 


Kt'J+1-<1J+'l<^  (18) 


for  any  u.r  <  j.    If  A3Jj  'J+    is  selected  in  any  given  path,  then  the  number  of  machines 
required  at  node  Vj+lj+1  is 

Pi+i  {dk+i  -  rUJ 


and  the  remnant  at  Vj+i,j+i  is 


A 


(19) 


rj+hj+i=  3-^^-dj+1-ruj.  (20) 

Pj+i 

Clearly,  the  cost  of  the  h-arc  Aji  will  vary  depending  upon  which  f-arc  was  selected  to 
reach  Vj+\j+\.  Thus,  the  single  arc  shown  in  Figure  1  actually  represents  a  set  of  arcs  with 
different  costs  although  these  costs  are  within  F2  of  each  other.  Clearly,  the  state  space  for 
the  shortest  path  problem  grows  exponentially  in  the  number  of  products. 

3.1.     A  Greedy  Algorithm 

If  remnant  differences  among  various  h-  and  f-arcs  incident  upon  the  same  node  are  ignored, 
then  P2  can  be  formulated  as  a  dynamic  program  in  which  the  stages  are  the  consecutively 
numbered  products,  and  the  states  at  a  given  stage  i  are  the  lowered  numbered  products  to 
which  product  i  can  be  assigned.  The  corresponding  network  is  identical  to  Figure  1. 

Define  g(t)  to  be  the  minimum  cost  of  reaching  stage  t,  t  =  1, . . . ,  N.  The  recursive  relation- 
ship is 

g(t)  =  minj<t{gjt},      where  (21) 

9jt     =    9j,t-\  +cj,'«-i'     J  <  l  (22' 

=    rnink<t  {gk,t-\  +  4!t-i}  ' 23  ' 

Hereafter,  we  refer  to  this  heuristic  solution  method  as  Greedy.  Clearly.  Greedy  requires 
0(N2)  computation  time. 
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Proposition  3.    Let  Z*  be  the  optimal  solution  value  to  P2,  and  let  ZG  denote  the 
solution  value  of  Greedy.   Then,  ZG  <  Z*  +  F2. 

PROOF:  Let  irG(-K*)  be  the  path  followed  by  Greedy  (the  optimal  solution).  Using  super- 
scripts LG'  and  '"'',  to  distinguish  the  variables  under  these  two  paths,  compare  the  values  of 
gG  and  g*t.  Note  that  the  two  solutions  can  differ  only  if  the  f-arc  selected  by  Greedy  to 
reach  node  Vjj  (say)  is  not  the  same  as  the  f-arc  selected  by  the  optimal  solution.  Hence, 
gG  =  g*t,  for  j  <  t  and  t  <  2.  Let  t  =  r  be  the  first  stage  where  these  values  differ.  Then 

gftT-i  =9j,T-v     Y;<T-i,     i<r-l 
and,  therefore,  by  induction 

»&_!  =  r*T_1?     Vj<r-i,     i<t-\  (24) 

From,  (30),  we  then  have 

<&  =  <£T,      V;<r  (25) 

It  follows  that  gG    =  g*T,  and  from  the  construction  of  Greedy, 

fr.r    <  g'r.r  (26) 

which  implies  that  ttg  and  7r*  use  different  arcs  to  reach  VTT  from  stage  r  —  1.  Let  TTG(ir*) 
use  arc  A^T_l(AlJT_l).  Consider  the  remnants  rUjT_i  and  rVjT_i  at  nodes  Vu,T-\  and  K,r-i- 
If  rU(T_!  >  rVfT_i,  then,  the  cost  of  any  arc  coming  from  VTT  in  7rG  can  be  no  more  than 
the  cost  of  the  corresponding  arc  in  7r*.   From  (26),  it  follows  that  <7?T+1  <  9j,r+v    Hence, 

gG(T+l)<g*(T  +  l). 

Now  suppose  that  rUi7._i  <  rw,T_i.  Of  the  two  arcs  coming  out  of  VTT,  consider  the  f-arc 
^t+i.t+i  £rst>  j£  ^g  remnant  at  VTT  is  ignored,  then  the  (possibly)  fractional  number  of 
machines  required  at  VT+itT+i  is 

_  pT+i"r+i  ,    a 

^T+l,r  +  l   —   "7 —  <^t  +  1,t  +  1  "T  Pr  +  l.r  +  l, 

where  aT+liT+1  and  /?r+i)T+i,  respectively,  denote  the  integer  and  the  fractional  part  of 
/iT+iiT+i.    Now  consider  the  impact  of  remnants  rUiT_i  and  rV(T_i  available  under  ttg  and 

12 


7r*  on  the  cost  of  arc  <4£*1,T+1.  Clearly. 

<*T+1,T+1     <     MJ+lj  +  1     <    AfV+lJ+l     <    CtT+l,T+l    +    1. 

If  -^h-i  ,+i  —  ^7+ij+n  then  the  remnant  differences  do  not  result  in  any  machine  savings, 
and  <??+1, r+1  <  g"r+i,r+v  From  (20),  it  follows  that  r«+lfT+1  -  r£+lfT+1  =  r»T  -  r*T.  This 
implies  that  the  remnant  differences  available  at  VTT  is  carried  forward  to  K.+itT+i. 

On  the  other  hand,  if  A/J'+1  +1  =  M'J  +  lj  +  l  +  1,  then,  remnant  differences  result  in  saving 
one  machine  at  V'r+i,r+i,  and 

r+l,r+l/      G\  „t+1,t+1'G\  T? 

Therefore,  gG+l<T+l  <  g'+lr+l  +  F2.  However,  in  this  case,  from  (20),  it  follows  that 


r11  -rv  -  — 

'  r+l,r+l  '  r+l,T+l   ~ 

Pr+\ 


+  P«    _  r"    >  0 


which  implies  that  xG  dominates  tt"  from  stage  r  +  2  onwards. 

In  summary,  either  i)  gG+l  T+l  <  j*+lT+1,  or  ii)  gG+1  T+1  <  5*4-1  t+i  +  ^21  and  ^  dominates 
t*  from  stage  r  +  2  onwards.  By  considering  the  h-arc  A£jT+1,  it  can  similarly  be  shown  that 
either  i)  gGT+i  <  9*,t+ii  or  n)  g?,r+i  <  #*.r+i  +  ^2,  and  7rG  dominates  7r*  from  stage  r  +  2 
onwards.  Clearly,  for  the  other  nodes  V^,T+i,  j  <  r  at  stage  r  -f  1,  we  have  gGT+i  =  9j,T+i- 
it  follows  that  either 

gG(r+i)    <    gm(r  +  l)t     or 
gG{r+l)    <    gm{r  +  l)  +  F2,     and 

gG(t)    <   g"(t)  +  F2,    t>r  +  2 

This  proves  the  proposition.  □ 

Corollary  2.    Let  Zt"   and  ZG  be  the  optimal  cost  and  the  cost  under  Greedy  to  reach 
node  VtJ  from  Vn.    Then,  Z"  >  Zfi  —  F2. 

3.2.     Exact  Solution  Algorithm 
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We  combine  corollary  3  with  Greedy  to  construct  an  improvement  algorithm  that  solves 
P2  exactly.  Given  the  Greedy  solution,  this  algorithm  starts  with  stage  N  and  construct  a 
search  tree  that  includes  selected  nodes  at  the  preceding  stages  in  G.  Hereafter,  we  use  the 
term  'vertex'  to  distinguish  a  node  generated  in  the  search  tree  from  a  node  in  G.  The  major 
aspects  of  the  algorithm  are  the  initial  solution  and  pre-processing,  and  branching,  updating 
and  fathoming. 

Algorithm  Exact 

Initial  Solution  and  Pre-Processing 

Solve  P2  using  Greedy.  Set  the  current  upper  bound  U B  equal  to  the  value  of  the  Greedy 
solution.  Determine  r^,  the  remnant  available  at  node  VlJ,  Vj,  and  Vz  <  j.  At  stages 
j  =  2,...,  iV,  compute 

Pj  (r*".i-i  _  rhj-\)  i 

Pij     =     -: ,    and 

Kij     =     [gij-i  +  cft-i)  "  {dh,j-i  +  4J,j-i) 

where  h  =  arg  mink<]  \gk,j-\  +  4fi-i}- 

Determine  Tj  =  {Vitj-i\Kij  <  F2)  for  j  =  2, . . . ,  T(=  N  -f  1).  Rank  nodes  within  each  Vj  in 
the  nonincreasing  order  of  kxj. 

For  each  node  Vjn  in  Tj,  compute  [Ijn  =  a3j^  +  /3jn. 

Branching,  Updating  and  Fathoming 

Construct  a  search  tree  rooted  at  T.  Generate  a  vertex  at  the  first  level  corresponding  to 
each  node  in  Tj-  For  each  such  vertex  v  that  corresponds  to  (say)  node  Vij  in  G,  compute 

@v  =  PjN\     and  Vv  =  9jn 

and  generate  vertices  at  the  second  level  corresponding  to  nodes  in  Tj.  Continue  generating 
nodes  in  this  fashion  until  either  a  node  is  fathomed  (to  be  described  below)  or  node  Vn  is 
reached.  For  any  vertex  v  at  level  2  or  below,  compute 

0v  =  0upv/pu,    and   zv  =  plJ+lpk/pJ+l 
14 


where  u  is  the  parent  vertex  of  u,  and  Kj(K',j+i)  is  the  node  in  G  that  corresponds  vertex  v 
(u). 

If  zv  >  #„,  then  fathom  vertex  r.  on  the  basis  of  remnant,  with  solution  value  tyv  —  >/n  ■+■ 
Kjj  —  F>.  If,  t'L,  <  UB,  then  set  c/i?  =  ipv,  and  record  u  as  the  incumbent.  Otherwise,  if 
zv  <  0,  ,  then  set  r,  =  i.-n  +  K{j.   If  tpv  >  t/B,  then  fathom  u  on  the  basis  of  cost. 

Let  v  be  the  incumbent  vertex,  which  corresponds  to  node  Kj,  at  the  end  of  the  procedure. 
Construct  the  optimal  solution  as  follows:  Trace  the  path  leading  from  v  to  the  root  node, 
and  find  the  nodes  in  G  corresponding  to  each  vertex  in  this  path.  These  nodes  determine 
the  corresponding  path  in  G  from  Vt]  to  T.  Find  the  shortest  path  from  V\\  to  V{j  using 
Greedy  to  complete  the  solution. 

Proposition  4.    Exact  solves  P2  optimally  in  0{N4)  time. 

PROOF:  Let  the  optimal  solution  have  pivots  (1),  (2), . . . ,  (n).  Then  from  corollary  3,  it 
follows  that 

Z"  <  ZG  <  g(n)N  <  Z*  +  F2. 

From  the  definition  of  Ij,  it  follows  that  Vtn\^  G  Tj.  Similarly,  it  can  be  shown  that 
V(n-\),(n)-\  G  r(n),  and  in  general,  V(n_,)(n_1+1)_1  £  F(n_1+1).  If  no  fathoming  takes  place,  all 
these  Ts  are  systematically  enumerated  in  the  search  tree;  clearly,  the  optimal  solution  must 
correspond  to  one  of  the  paths  in  this  tree. 

It  only  remains  to  show  that  fathoming  does  not  cut  off  the  optimal  solution.  Clearly,  cost 
based  fathoming  does  not  do  so.  In  remnant-  based  fathoming,  a  vertex  v  is  fathomed  if 
whenever  an  f-arc  leading  to  that  vertex  is  found  with  a  remnant  large  enough  to  absorb 
the  dangling  fractional  machine  on  the  last  line.  Clearly,  no  further  benefit  is  to  be  gained 
by  considering  vertices  that  are  children  of  v  because  no  more  machine  savings  are  possible, 
and  from  the  construction  of  Greedy,  the  cost  of  reaching  the  node  corresponding  to  v 
under  Greedy  is  the  minimal  among  all  paths.  This  completes  the  proof  of  the  optimally 
of  Exact. 
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The  number  of  vertices  at  the  first  level  for  —  \^t\  <  N.  A  given  vertex  V3n,J  >  2  at 
this  level  will  generate  \Tj\  <  j  —  1  children,  and  so  on.  In  the  worst  case,  the  number  of 
descendants  generated  by  Vjjv  is 

4>jn    =    [(j-l)  +  (i-2)  +  ...  +  2]  +  [(i-2)  +  (i-3)  +  ...  +  2]  +  ...  +  [2] 

q-i)0"-2)(j  +  +3) 

6 
Since  there  are  N  such  vertices  at  level  1,  the  overall  number  of  vertices  generated  is  0(N4). 
At  each  vertex,  there  are  2  computations  and  1  comparisons  done.  Hence,  the  overall  com- 
plexity of  the  algorithm  is  0(N4).  □ 


4      No  Splitting 

The  unused  capacity  at  line  i,  termed  the  remnant  at  line  i,  is  given  by 

r*  =  —  -  E  *■ 

/".    <«, 

Ri  is  the  additional  number  of  units  of  any  product  j,j  >  i  that  could  be  processed  at  line 
i.  Ideally.  R(  should  be  zero;  however,  because  the  demand  of  any  product  cannot  be  split 
across  more  than  one  line,  Ri  is  likely  to  be  positive  in  general. 

Consider  a  subgraph  S  C  G  that  represents  assignments  of  products  1  through  A',  K  <  P. 
Let  7r  be  any  path  in  S  with  pivots  /j, . . . ,  /n,  and  let  £,-,  i  =  1, . . . ,  n  be  the  products  assigned 
to  line  i.  Consider  an  identical  product  assignment  represented  by  path  tt'  in  G' .  In  tt' ,  the 
number  of  machines  required  at  line  i,i  >  2  is  given  by 

Pi,  (Utec,  dt  ~  Ri-i) 


n'  = 


A 


=     rii 

<     nt.  (27) 


Hence,  demand  splitting  can  lead  to  possible  savings  in  the  number  of  machines  required. 
Let  the  total  number  of  machines  saved  in  tt'  over  tt  be  a.  The  capacity  available  at  line 
i  for  possibly  absorbing  demands  of  products  K  +  1  through  P  is  Rt  =  i?,-,  and  the  total 
capacity  available  is 

n  n  n       a      .  n 

t'=l  i  =  l  i=l    rU  i=l  tec, 
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In  k\  the  capacity  R't  available  at  line  i,  i  ^  n  equals  zero,  and 

An,, 


K  =  K   = 


vu 


n-l 


tecn  1=1 


=   tn~-  (28) 

,=  1  P'n 

Let  Nuv  denote  the  number  of  machines  required  by  product  v  on  line  u.  Then  Nuv  =  \7uv]i 

whore 

Pu{dv  -  ru) 

and  ru  is  the  remnant  available  on  line  u  based  upon  the  set  of  products  Cu  already  assigned 
to  it.  Clearly,  if  u  =  v,  then  Cu  =  0,  and  ru  =  0.  The  cost  of  arc  Af?  is  then  given  by 

c™    =    Fl  +  F2NUV     \tu  =  v 
—     F2NUV     otherwise 

Note  that  the  costs  of  all  arcs  incident  upon  Vuu  are  the  same.  However,  the  costs  of  arcs 
incident  upon  Vuv,  u  ^  u,  may  vary  because  they  depend  upon  ru  which  in  turn  depends 
upon  Cu.  Let  juv  —  Nuv  +  fuv,  where  Nuv  =  [luv\  and  fuv  is  a  fraction.  We  then  have  the 
following  result: 

Proposition  5.  Let  it'  be  an  optimal  path  from  Vu  to  t  in  G.  If  a  b-arc  ,4"J  lies  on  this 
path,  tli (ii  A'IJ  +  i  =  Nu<J+i. 

4.1.     Sequential  Assignment  Policy 

liilike  the  demand  splitting  case,  the  optimal  solution  to  PI  does  not  possess  the  sequential 
assignement  property.  However,  this  property  can  be  enforced  to  obtain  a  useful  heuristic 
solution.  Consider  the  following  problem: 

Tree  Partitioning  Problem  (Corneujols,  Nemhauser  and  Wolsey  1990):  Suppose  G  —  (A  .A) 
is  a  tree  graph,  and  A  is  a  \Af\  by  \J\f\  matrix  with  elements  Dl3  that  represent  the  distance 
betweeD  nodes  v,  and  r,,  for  all  u,-,  v}  G  A*.  Let  the  weight  of  any  subtree  Gj  =  (A  j.A,\ 
be  w(@j)  =  maxv  --ft  [HveAf  ^ti)- Determine  the  partition  of  G  into  subtrees  such  thai 
minimizes  the  sum  of  subtree  weights. 
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Corneujols  et  al.  (1990)  present  a  dynamic  programming  algorithm  that  solves  the  tree 
partitioning  problem  in  OdA/'l2)  time. 

Given  this  property,  the  sequential  assignment  policy  can  be  modeled  as  a  tree  partitioning 
problem  as  follows.  Let  Af  =  {v\,. . .  ,vn}  denote  the  products  1  through  N,  and  let  A  = 
{Ay.+l  :  i  =  1, 2, . . .  ,N  —  1}.  There  is  a  unique  path,  comprising  one  or  more  pivots,  from 
v i  to  uyy.  Furthermore,  in  any  solution  to  this  problem,  each  pivot  generates  a  subtree.  The 
equivalence  is  complete  if  we  assign 

diPi 


Da  =  F1  + 


A 


and 


Da     =     F, 


u 


Pi  EJ=j  dt 


A 


if  i  <  j 


=    oo      otherwise 


It  follows  that  the  optimal  SAP  solution  can  be  found  in  0(N2)  time.  It  is  useful  to  note 
that,  under  SAP,  the  cost  of  all  arcs  incident  on  node  Vuu  is  the  same.  In  addition,  because 
SAP  considers  only  the  f-arcs  and  b-arcs,  we  have 

Remark  1.    If  there  are  no  b-arcs  in  any  optimal  solution  to  PI,  then  it  is  given  by  the 
sequential  assignment  policy. 

5      Lower  Bounds  on  PI 

We  propose  four  lower  bounds  which  are  discussed  below. 

5.1.     Fractional  Machines 

The  first  bound  LB1  is  obtained  by  relaxing  constraint  (6)  to  read  n3  >  0.    In  that  case, 
constraint  (3)  will  be  satisfied  as  an  equality.  In  particular, 

rij  =  ^JcjjXj'j,      where  C{3  —  diPj/A. 


PI  then  reduces  to 
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P1F 

Z  -  min  ^  Fiyj  +  ^  ^c^jx^j 
j  j 

subject  to 

Z  -ru  =  1 
i€Ji 

Xij  <  Vj 
Xij,Vj  €  {0.  1}. 
Proposition  6.     There  exist*  an  optimal  solution  to  PlF  that  follows  SAP. 

PROOF:    We  show  that  in  at  least  one  optimal  solution  to  PlF.  if  product   /.  /   £  A    is 

produced  on  line  j.j  <  /.  then  j  +  1 /  —  1  are  also  produced  on  j.  Let  o~\  be  an  optimal 

solution  to  PlF  that  does  not  have  this  property.  We  show  that  an  alternative  solution  o~2. 
no  worse  than  c^,  that  does  have  this  property. 

From  our  assumption,  there  has  to  be  at  least  one  product  t,j <  <  t  <  i '  —  1  that  is  produced 
on  line  fc,  k  ^  j  in  a\.  Clearly,  k  <  t  <  i.  Let  z(-)  denote  the  cost  of  solution  (•).  Consider 
the  following  cases: 

Case  I:  k  <  j 

Let  o~2  be  identical  to  o~\  except  in  that  t  is  produced  on  line  j.  Then 

i      \         t      \       dtPk       dtpJ        (  \dt  ^  n 

Z(<?l)  -  Z(<T2)  =  -^ £"  =  (Pk  -Pj)-£    ^°- 

Case  II:  k  <  j 

In  this  case,  let  a2  be  identical  to  o~\  except  in  that  products  t  + 1, . . .  ,  i  are  assigned  to  pivot 

k.  Then 

P;(ZUirf7)       P*  (EU+i  ^) 
Z((7!)  -  2(a2)  = >  U. 

Repeating  this  argument  for  all  such  t  gives  the  desired  result.  □ 

LB1  can,  therefore,  be  determined  by  an  0{P2)  dynamic  programming  algorithm  as  discussed 
in  §2. 

5.2.     Demand  Splitting 
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The  second  bound  is  obtained  by  using  a  greedy  algorithm  for  solving  the  demand  splitting 
problem.  Let  the  solution  obtained  by  this  approach  be  LB2. 

Proposition  7.    LBS  is  a  valid  lower  bound  on  the  optimal  solution  value  for  Pi. 

PROOF:  We  show  that  the  cost  of  the  optimal  path  ir*  in  G  from  V\\  to  T  is  no  less  than  the 
cost  of  at  least  one  path  between  these  two  nodes  in  G'  that  is  considered  under  Greedy. 

Let  a  be  a  segment  of  7r*  between  nodes  Vu  and  VX]  such  that  it  consists  only  of  f-  and 
h-arcs.  Note  that  such  a  segment  does  exist  since  the  arcs  emanating  from  Vu  consist  of  one 
f-arc  and  one  h-arc.  Also  note  that,  because  it  does  not  include  any  b-arcs,  a  exists  in  G'  as 
well  and  is,  therefore,  considered  by  Greedy.  Let  r/(-)[r]'{-)}  denote  the  number  of  machines 
required  by  segment  (•)  in  G{G'   under  Greedy).  Then 

Remark  2.    n(a)  —  n'(o-)  =  a  >  0. 

Noting  that  the  number  of  pivots  in  a  is  the  same  in  both  G  and  G'',  and  taking  Vij  =  T  in 
Lemma  1,  completes  the  proof  of  the  proposition  if  ir*  consists  only  of  f-and  h-  arcs.  Next 
consider  the  case  in  which  ir*  consists  of  one  or  more  b-arcs.  Let  the  node  at  the  tail  of  the 
first  b-arc  be  Vuv.  Let  a\  denote  the  first  segment  on  tt*  that  starts  with  this  b-arc  and  ends 
with  an  f-arc  incident  upon  a  node  that  belongs  to  pivot  9  such  that  9  >  u.  Note  that  there 
will  always  be  such  a  segment  given  that  the  terminal  node  T  is  reachable  only  with  an  f-arc. 
Clearly,  o~i  has  to  be  preceded  by  at  least  one  f-arc,  and  therefore,  u  >  2.  Let  the  set  of 
pivots  visited  by  <7i  be  C2  =  {a/},  and  let  A\  denote  the  set  of  products  assigned  to  pivot 
/,  /  £  C2  prior  to  o~\.  G\  can  then  be  expressed  as  A^"1"1  —  •  •  •  —  Aa'v*+£    for  some  q  £  C2. 

Let  a0  denote  the  segment  of  7r*  that  starts  with  Vu  and  ends  with  the  arc  incident  upon 
Vuv,  and  let  the  set  of  pivots  visited  by  a0  be  C1.  At  the  end  of  a0,  let  n(cr0)  denote  the 
remnant  on  line  /,  /  £  Cl  in  G,  and  R'u(a0)  be  the  remnant  on  line  u  in  G' .  From  Corollary 
1,  C2  C  C1.  Also  note  that 

Vu  <  Pi,      V/  £  C2.  (29) 

From  remark  2,  a  =  n(cr0)  —  n'(ao)  >  0.  Now  consider  the  following  cases: 

Case  I:  9  —  u. 
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In  this  case,  we  show  that  rrx  is  dominated  by  the  segment  a2  in  G\  where  (j2 


nu.v 


Au,v  +  k       ■     N©W 


and 


rj'(a2)  = 


Pv 


\ 


7(*i)    =     E 

/€£2 

>    £ 

/€£2 


Pi 


> 


E/<E£2  P«  (E«M,  r/<  ~  r'j 


Pu  [(E/g£2  E<€.4,  ^)  -  E/g£2  n 


> 


Pu 


E/e£2  E<eA  ^)  ~~  ^u 


,\ 


—  a 


=     r)'(a2)-a 


Hence,  rj(cro)  +  t]{g\)  >  r/'(<7o)  +  T]'(a2).  Note  that  no  new  pivots  are  opened  under  both  (T\ 
and  rr2. 

Case  II:  6  >  u. 

In  this  case,  note  that  9  =  v  +  k  +  1.   We  show  that  cri  is  dominated  by  the  segment  a2  in 

For  the  first  k  arcs,  the  arguments  used  in  case  I  apply  here  as  well.   For  the  (p  +  l)st  arc. 
let  f)'{{n)  denote  the  number  of  machines  required  under  a2  in  G'  {(j\  in  G).  Then 


K  = 


Pejde-K) 
A 


< 


Pod  6 
A 


=  A 


Hence,  as  before,  //(cr0)  +  ^(^i)  >  v'i^o)  +  rfi^)- 

Put  q  =  77 ((To)  +  r)((j\)  —  r]'((T0)  +  f]'(cr2)-  Then,  we  have 


n.  1 

Pe 


21 


The  desired  result,  folio ws  by  repeating  the  above  arguments  for  each  such  segment  0\  in  G. 

a. 

The  following  result  shows  that  this  bound  is  stronger  than  LB\. 
Proposition  8.    LB2  >  LBl. 

PROOF:We  show  that  the  cost  of  any  path  tt  from  Vn  to  T  under  Greedy  is  at  least  as  large 
as  the  cost  of  the  same  path  when  fractional  machines  are  permitted.  Let  the  set  of  pivots 
on  7r  be  £.  and  let  L  =  \C\.  Also,  let  £/,/  6  C  denote  the  set  of  products  assigned  to  pivot 
/.  Then  the  total  cost  of  tt  given  fractional  machines  is 


C1(tt)  =  LF1  +  F2  £ 


PiL 


tec,  at 


l€calL 


A 


and  the  cost  of  it  under  Greedy  is 

C2(;r)    =    LFl  +  F2Y, 


/€£ 


=    LFl  +  F2Y, 


tec 


=    dW  +  l/A^-W 
lec 


=  CiW  +  iM 
>  cm*) 


H   R'iiPi  ~P/+i)  +  R'l 
liec\L 


(30) 


5.3.     Lagrangean  Relaxation 

The  third  bound  LBS  is  obtained  by  relaxing  the  problem  by  associating  nonnegative 
multipliers  \j  with  constraints  (2).  The  resulting  problem  is 

PI  -  LR 

p  jP_       /  _  \ 

(31) 


p  p       l 

Minimize       ^(Fiyj  +  F2rij)  -  J2  ^     J2  x*i  ~  1 

j=i  i=i       \jeJ, 


subject  to 


(3) -(6) 
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For  given   multipliers.  Pi    -    LR  separates  into  P  independent  problems,   where  the  jth 

problem  is 


PI  -  LR, 


subject  to 


Minimize       F\\jj  +  F^rij  —  ^  \{X{j  (32) 


Pj  (£**«]  <Anh  (33) 

n3  <  BVj,  (34) 

ay  €{0,1},     zelj  (35) 

yj  6  {0,1};^  >  0.    integer,  (36) 

This  problem  can  be  solved  by  solving  the  following  knapsack  problem 

Minimize      Z  =  F\  +  F2TIJ  —  ^  A,x,_, 

subject  to 

Pj     Y,  dixH  I  ^  Ann 

Xij  e  {0, 1},    1  €  ij 

If  the  optimal  value  of  Z  =  Z*  <  0,  then  i/j  =  1,  else,  it  is  zero.  The  above  problem  can  be 
converted  into  a  standard  knapsack  formulation  by  substituting  x'tJ  =  1  —  xl3  to  yield  the 
formulation 

Minimize      Z  —  F\  +  F2rij +  ^  ^ix'tJ  ~  ^2  ^* 

subject  to 

Pj  I  Z  fi'xu  )  ^  ft     Z  rf«  ~  AnJ  1  - 
x,j  E  {0,1}.     1  £lj 
The  subgradient  optimization  technique  is  used  to  adjusting  the  multipliers. 
5.4.     Surrogate  Relaxation 
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The  fourth  bound  LBA  is  obtained  by  associating  nonnegative  multipliers  X3  with  constraints 
(3)  and  forming  a  surrogate  constraint.  The  resulting  problem  is 


PI  -  SR 


subject  to 


p 
Minimize       22(FiVj  +  F2rij)  (37) 


^Xij  =  l,    ietf  (38) 

J2  Yl  XjPjdixiJ  <  H  XjArij,    j  e  Af  (39) 

nj<Byh    jetf  (40) 


x 


'j 


e{o,i},    ijeAf  (41) 


yj  e  {0, 1};  rij  >  0,    integer,     j  G  .V  (42) 

The  surrogate  constraint  essentially  insures  that  the  slack  capacity  available  at  any  line 
Sj  =  An  j  —  pjdtXij  is  scaled  by  the  multiplier  \31  and  is  added  to  the  capacity  of  the  next 
line.  Clearly,  if  A^  >  1,  then  the  capacity  added  is  greater  than  the  slack  available.  Hence,  in 
order  to  insure  that  a  tight  bound  is  obtained,  we  restrict  ourselves  in  the  multiplier  space 
to  \j  <  1.  Also,  note  that  if  \3  =  pJ+i/pj,  then  the  amount  of  slack  capacity  that  is  passed 
from  line  j  to  line  j  +  1  is  capable  of  producing  exactly  r2  units,  and  therefore,  the  optimal 
solution  to  PI  -  SR  is  the  same  as  the  optimal  solution  to  the  demand  splitting  problem 
P2.  Because  we  know  that  the  Greedy  solution  that  is  an  upper  bound  to  P2  is  a  lower 
bound  to  PI,  our  objective  in  this  relaxation  is  to  generate  multipliers  that  result  in  tighter 
bounds  still.  However,  the  above  values  of  these  multipliers  can  serve  as  good  intial  values. 

Proposition  9.     The  sequential  assignment  policy  is  optimal  for  P1A. 

PROOF:  The  proof  is  similar  to  that  of  proposition  3,  and  is,  therefore,  omitted. 
Multiplier  adjustment  is  based  on  the  subgradient  optimization  method. 
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6  An  Implicit  Enumeration  Algorithm 

An  exact  solution  algorithm  for  solving  Pi  that  is  based  on  implicit  enumeration  can  now  be 
designed.  The  root  node  in  the  search  tree  corresponds  to  product  1,  and  each  subsequent 
level  corresponds  to  successive  products  considered  in  the  order  of  their  indexes.  At  a  given 
level  /,  the  leftmost  branch  represents  the  SAP  assignment  of  product  t.  In  addition,  a 
branch  is  created  for  every  b-arc  that  satisfies  proposition  5.  The  search  tree  generated  in 
this  manner  is  exhaustive  because  it  enumerates  each  undominated  b-arc.  Suppose  that 
a  branch  corresponding  to  b-arc  A3t'tt+  is  being  considered  for  further  augmentation.  The 
corresponding  graph  is  then  regenerated  by  combining  products  j  and  t  +  1  into  a  single 
product  with  processing  time  p3,  and  demand  dj  +  c^+i-  In  so  doing,  we  insure  that  the 
resulting  graph  contains  only  f-arcs  and  h-arcs.  The  leaf  nodes  then  correspond  to  different 
completions  of  graph  Q ,  corresponding  to  different  selections  of  eligible  b-arcs.  However, 
since  these  completions  contain  only  f-arcs  and  h-arcs,  they  are  optimally  solved  optimally 
solved  in  0(|A'|2)  time. 

The  lower  bound  LB  at  any  node  is  given  by 

LB  =  max  {LB2,  LBZ,  LB4) . 

However,  it  may  be  computation-effective  to  use  only  one  bound,  say  LB2.  The  upper  bound 
at  any  node  can  be  obtained  by  applying  the  sequential  assignment  policy  to  the  remaining 
subtree. 

7  Conclusion 

In  this  paper,  we  have  formulated  and  analyzed  the  problem  of  providing  several  manufac- 
turing lines,  and  assigning  products  to  these  lines.  We  consider  two  versions  of  this  problem. 
rhe  first  version  permits  splitting  the  demand  of  any  product  across  multiple  lines.  We 
analyze  the  structure  of  this  problem  and  develop  a  polynomial  time  algorithm  for  solving 
it.  Next,  we  consider  the  case  in  which  each  product  is  assigned  to  exactly  one  line.  We 
show  thai   this  problem  is  NP-hard.    We  explore  the  properties  of  the  optimal  solution  to 
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this  problem  as  well  as  the  Lagrangean  and  surrogate  relaxations.    Finally,  we  present  an 
implicit  enumeration  method  for  solving  this  problem. 
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